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ON  THE  C02  LASER  GENERATED  CHANNEL  CONDUCTIVITY  AND 
ELECTRON  BEAM  IONIZATION  OF  AMMONIA 


I .  INTRODUCTION 

Gaseous  channels  with  appropriate  conductivity  can  fascilitate  the 
transport  and  guidance  of  charged  particle  beams  and  electrical  discharges. 
Such  channels  can  be  generated  by  focusing  a  beam  of  high  power  laser  or 
microwave  radiation  to  breakdown  the  gas  over  the  desired  length.  The  gas 
breakdown  results  in  the  formation  of  a  plasma  and  the  heating  of  the  gas, 
thereby  generating  a  conducting  channel.  Studies1-3  have  often  used  C02  laser 
to  form  channels  in  various  gases  such  as  air,  N2,  NH3  and  C2H2.  The 
characteristics  of  such  channels,  i.e. ,  the  channel  temperature  and  the 
channel  conductivity  must  be  understood  in  order  to  provide  the  basis  for  the 
understanding  of  charged  particle  beam  transport  and  guidance  in  the  channel. 

In  this  report  we  present  a  preliminary  analysis  on  the  characteristics 
of  NH3  channels  formed  by  using  a  CO.,  laser  where  the  laser  power  is  below  the 
threshold  for  breakdown.  The  report  also  presents  the  electron  beam 
ionization  (collisional  and  avalanche)  of  NH3  and  calculates  the  conductivity 
generated  by  the  beam  for  two  typical  experiments  (NRL  and  Sandia 
Experiments).  The  calculation  of  the  electron  beam  generated  conductivity 
shows  that  the  beam  does  its  own  thing  in  that  the  residual  channel 
conductivity  plays  no  role  but  a  neutral  density  channel  left  behind  by  the 
laser  will  modify  the  conductivity  generated  by  the  beam.  This  is  because  the 
channel  conductivity  is  either  non  existent  or  below  the  relevant  value  in 
these  two  experiments. 

II.  C02  LASER  ABSORPTION  BY  NH3< 

The  C02  laser  is  a  multi-line  laser,  and  in  general  it  can  be  tuned  to  a 
given  frequency.  Several  of  the  laser  lines  are  in  near  resonance  with 
absorption  bands  in  NH3»  As  an  example,  the  \>2  mode  energy  diagram  is  shown 
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in  Figure  1  where  the  C02  laser  transitions  i.e.  R(16)  and  R(30)  are  shown  to 


coincide  with  transitions  within  the  v*2  band  of  NH^.  Such  punping  in  general 

4 

also  leads  to  lasing  within  NH^  at  frequencies  smaller  conpared  to  the  pump 
frequency.  Hie  absorption  of  the  laser  quanta  by  the  molecule  in  the  ground 
state  would  raise  the  molecule  to  V=1  level  and  several  other  processes  of 
inportance  will  occur. 

1.  The  vibrational  energy  is  converted  to  gas  kinetic  energy  through 
the  VT  energy  exchange. 

2.  Higher  levels  can  be  excited  by  further  absorption  of  the 
vibrational  quanta,  resulting  in  collisionless  dissociation  of  the 
molecule. 

3.  Vibrational-vibrational  collisions  resulting  in  the  ladder  type 
excitation  and  dissociation. 

4.  Stimulated  emission. 

5.  Saturation  of  absorption. 

C 

Hie  first  process  is  very  rapid  and  has  a  relaxation  time  of  2.5  nsec/atm. 
Hie  second  process  occurs  for  short  pulse  lasers  ~  in  the  n  sec  regime,  while 
the  last  process  occurs  for  lasers  of  long  duration  ~>  10  nsec. 


t§ 


III.  THE  NRL  EXPERIMENT 

The  experiment  at  NRL  carried  out  by  the  Channel  Physics  Group  in  the 
Plasma  Physics  Division  utilizes  a  home  made  maltimode  C02  laser.  This  laser 
has®  a  total  power  of  (50J)  and  pulse  length  of  200  nsec.  It  was  used  to 
study  channel  physics  in  NH.j  where  a  40  Ttorr  NH^  filled  the  chamber  through 
which  the  C02  laser  was  passed  in  weak  focus  condition  (The  focusing  was 
outside  the  chamber).  Hie  energy  absorbed  by  the  gas  was  25  mJ  per  cm^ 
resulting^  in  a  tenperature  rise  to  ~  0.1  eV. 
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IV.  THERMAL  DISSOCIATION  AND  IONIZATION  OF  NH3 

A.  In  this  section  we  will  discuss  the  dissociation  and  ionization  of 
NEj  based  on  the  equilibrium  calculation.  The  absorption  of  C02  laser 
radiation  by  NH3  and  the  rapid  relaxation  of  the  vibrational  energy  to  kinetic 
energy  justifies  the  equilibrium  approach.  This  relaxation  time  t  has  been 
measured^  in  pure  NH^  and  was  found  to  be  equal  to  2.5  nsec/atm.  This  time  in 
the  case  of  the  NRL  experiment  is  47.5  nsec  which  is  still  much  smaller  then 
the  pulse  width  (  ~  200  nsec.). 

Q 

The  mass  action  law  for  the  dissociation  equilibrium  is 


H  NH2 


Where  K  is  the  equilibrium  constant  H,  NH2  and  NH^  are  the  densities  of  the 

Q 

relevent  species.  The  equilibrium  constant  can  be  expressed  as 


Q(H)  Q(NH2) 
Q< nh3) 


- ae/rt 


Where  Q(S)  is  the  total  partition  functions  of  species  S  and  is  a  function  of 

Q 

tenperature  only.  The  total  partition  function  in  general  is  a  product  of 
translational  and  internal  partition  functions  of  species  3  i.e. ,  Q  =  Q. 


We  shall  assume  that  the  major  contribution  to  the  internal  partition 
functions  of  the  molecule  is  due  to  the  rotational  and  vibrational  excitation 


which  is  always  the  case  for  poly-atomic  molecules  at  tenperatures  below 
lev.  The  rotational  partition  function  for  NH2  and  are8 


Q  =  0.00693  x 
r 


10 


V 


3 

T  (I 


aW 


(4) 


Where  IA,  IB  and  Ic  are  the  moments  of  inertia  along  the  principal  axis  or 


planes.  On  the  other  hand,  the  vibrational  partition  function  is 

.1 
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«v  - 


di 
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Where  ui^  is  the  frequency  of  the  vibrational  excitation  and  d^  is  the  degree 
of  degeneracy.  Utilizing  these  relations  into  the  equilibrium  constant  one 
obtains 


V 


d. 

1 


„  -w.hc/kT, 

3/2  21T  nh2  *  1  NH2 

K  =  ( 0 .3  IT)  x  10 


(Wc)„_  _  -<u.  he  AT 

NH3  V1"  1  )NH3 
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The  vibrational  partition  function  is  calculated  for  several  temperatures  and 
are  shown  in  Table  I 

TABLE  I 


T(°K) 

1000 

2000 


2V(NH2) 
1.26 
1 .92 


2.07 

5.9 


The  results  in  Table  I  are  obtained  using  the  following  excitation  energies' 


0.4,  0.12,  0.43  and  0.2  for  the  vibrational  levels  of  UH,  and  0.4,  0.19  and 


0.4  for  NHj.  Assuming  that  the  moments  of  inertia  of  and  NH-j  are  the 

same,  we  obtain 


3/2  21 

K  =  (0.3T)  /  x  10 


(7) 


Where  R  is  the  ratio  of  the  vibrational  partition  functions  of  NH2  to  NHj.  If 
we  designate  a  as  the  degree  of  dissociation  i.e.  a  =  then  the  mass  action 
law  can  be  written  as 


2 

a  MH3  =  K 


3 /2  21 

(0.3T)  x  10 


R( T)  Exp 


/  4983p\ 

V  / 


(8) 


Values  of  a  for  several  temperatures  are  given  in  Table  II  for  NH-j  =  1.4  x 
1018  cm"3  (  ~  40  Torr) . 


TABLE  II 


T(  °K) 

a 

1160 

1 . 2  ( —  6 ) 

2000 

2(  -2) 

2200 

1(-1) 

B.  The  thermal  ionization  of  in  most  probability  arises  as  a  result  of 

the  associative  ionization.  i.e. 


H  +  NH2  +  NH3  +  e 


(9) 


rather  than 


nh3  +  nh3  *  nh3  +  mh3  +  e 


(10) 


We  shall  calculate  both  rates,  however,  to  delineate  the  most  favorable 
regime. 

Saha's  Equation  for  the  last  process  is 


2  3  ,2 

Ne  =  2  (2jrmkT)  (- 1 17740/kT) 

NH3 


(11a) 


Ne  „  „  ,„15  „3/2  f  1  17740\ 

—  =  4.8  x  10  T  Exp  (- 


(11b) 


where  we  have  assumed  equal  statistical  weights  for  NH3  and  NH3+.  The 
equilibrium  electron  density  is  given  in  Table  III  for  several  tenperatures  of 
interest  for  an  NH^  density  of  40  Torr. 


TABLE  III 


T(  °k) 

1 160 
2000 
2200 


_3 

N  (cm  ) 
e 


1.5  x  10 


1.6  x  10c 


2.8  x  10' 


Mow  let  us  consider  the  associative  Ionization  process.  For  this  we  shall 


utilize  the  general  cross  section  developed  by  Melson  and  Dahler  . 
threshold  expression  for  this  cross  section  is 


The 


ax 


2  •  1  R  X  1  3 

*  -J 

— -  (E  -  Eth)  ’ 


(12) 
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Where  Rc  is 

the  crossing  point. 

Eth 

is 

the 

threshold 

for  association  and 

v, 

X  -  2(2u)  2 

. 

T  /E,  .  Here,  f  is 

the 

width  of 

. 

the  transition  and  E_  is  the 

c  d  c 

d 

derivative 

of  the  potential  c 

urve 

at 

the 

crossing 

point.  Using1  ' 

-3  - 
r  ^  0.04  and  R  ~  10  cm  i.e.  x  ~  1*6  x  10  and  E_  ~  1.0  x  10  ,  we  obtain 

c  ~  c  —  d  - 


a  . 

ai 


4.2  x  10 


16 


3/2 

(E  -  5.8) 

E 


(13) 


where  E  is  in  units  of  eV.  Using  the  slope  of  this  cross  se<  on  one  obtains 
the  following  rate  coefficient  for  the  associative  ionizatior 

-12  _5  8 /T 

R  .  =  10  /T  (5.8  +  2  T  )  e  *  '  g  (14) 

ai  g  g 

where  Tg  is  in  units  of  eV. 

If  we  assume  that  the  dissociative  recombination  rate  coefficient  is 
—  6  3 

10  cm  sec,  then  at  equilibrium  we  will  have 

-6  2 

10  N  =  R  .  H  NH2  ( 15) 

6  ax 


Table  IV  gives  the  values  of  the  electron  density  using  the  equilibrium  values 
of  H  given  in  Table  II 


( * )  Values  relevent  to  N  +  O  +  NO  +  e 


TABLE  IV 
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N  (cm  ) 

T(  °K) 

H  (cm-3) 

e 

1  160 

1.68  x  1012 

5.3  x  10~4 

2000 

2.8  x  1016 

2x10° 

2200 

1.4  x  1017 

5  x  107 

If  one  lowers  the  dissociative  recombination  rate  coefficient  by  an  order  of 
magnitude,  the  equilibrium  electron  density  will  be  as  in  Table  V. 


TABLE  V 


_3 


N  (cm  ) 

T(  °  K) 

6 

1150 

1.8  x  10“3 

2000 

6.3  x  106 

2200 

1.5  x  10® 

V.  THE  AVALANCHE  IONIZATION  OF  NH-j. 

The  Townsend  ionization  coefficient  for  NH^  has  been  measured  by  Risbud 

and  Naidu10  for  —  between  115  x  10~17v-cm2  and  2060  x  10-17V-cm2.  Tie 
N 

following  expressions  fit  the  data  to  better  than  30%. 


a 

N 


1.07  x  10 


Exp 


1  .23  x  10 
E/N 


) 


(16) 


for 


1 15(  -  17) <E/N  <  182  (-17) 


and 


(17) 


a  . -lb  _  ,  65.2  x  10~  . 

N  =  4*34  *  10  ^  ( - m ’ 


1 82( - 17)  <  -  <  2060  (-17) 


The  electron  drfit  velocity,  however,  is  measured1 1  over  a  limited  range,  i.e. 

—  1  7  2 

up  to  E/N  *  50  x  10  V-cn  .  We  have  obtained  the  following  fit  for  the 


7  »  -  1.26  x  10  t  2.0  x  10  (E/N) 

a 


drift  velocity  for  the  region  of  —  from  46.1  x  10-17  7-cm2  to  54.2  x  10~17  V- 

N 

2 

cm  .  This  expression  can  be  extended  to  higher  E/N  until  further  data  becomes 
available. 

Using  Equations  16  -  18  we  obtain  the  following  expressions  for  fthe 


reduced  ionization  frequency  in  NH-j 


,22  , E. 


—  =  [-1.26  x  10  +  2  x  10  ( — )  ]  x  2.0  x  10  Exp  (■ 

N  N 


1.23  x  10 


-)  (19) 


115  (-17)  <  E/N  <  182  (-17) 


S  22 

[-1.26  x  10  +  2  x  10  (E/N)]  4.34  x  lo'  Exp  (-  °5 1 0 - ) 


182  (-17)  <  E/N  <  2060  (-17) 


The  reduced  ionization  frequency  ^  is  given  in  Table  7I« 
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TABLE  VI 


E/N 


115  (-17)* 
123  (-17) 
154  (-17) 
200  (-17) 
216  (-17) 
308  (-17) 
462  (-17) 
617  (-17) 
925  (-17) 


0.494  (-1  1) 
1.14  (-11) 
1.06  (-10)  . 
6.50  (-10) 
8.85  (-10) 
3.08  (-9) 
9.56  (-9) 
1.82  (-8) 

3.9  (-8) 


1234  (-17)  6.17  (-8) 

1543  (-17)  8.64  (-8) 


VI.  THE  ELECTRON  BEAM  COLLISIONAL  IONIZATION  OP  NH3 

The  electron  beam  ionization  of  NH^  can  be  estimated  readily,  provided 

one  knows  the  NH3  stopping  power  for  higfr  energy  electrons.  For  this  purpose 

we  shall  assume  that  ammonia  is  neonlike  and  its  density  is  equivalent  to  that 

1  2 

of  nitrogen  atom.  Using  the  calculated  data  for  dE,/dx  for  neon  one  obtains 
a  value  of  1.96  x  10  eV/cm,  in  one  atmosphere  of  NH^,  for  electrons  with 
energy  of  1  MeV.  Thi3  value  increases  by  ~  20%  when  the  beam  electron  energy 
is  10  MeV. 


^ ^  (-17)  indicates  10-1^ 
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VII.  CALCULATIONS  AND  CONCLUSIONS 

In  the  preceeding  sections  we  have  presented  the  electron  beam 
collisional  and  avalanche  ionization  parameters  of  NH^.  Furthermore,  we 
presented  the  equilibrium  calculations  of  an  NH-j  channel  generated  by  lasers 
in  the  teaperature  range  of  0.1  to  0.19  eV.  From  these  equilibrium  relations 
we  have  obtained  the  residual  electron  density  in  the  same  tenperature  range. 

In  this  section  we  would  like  to  calculate  the  residual  conductivity  for 
two  channels,  one  at  a  temperature  of  0.1  eV  (NRL  Experiment)  and  the  other  at 
a  tenperature  of  0.19  eV  (Sandia  Experiment).  Then  we  will  calculate  the 
conductivity  generated  by  two  electron  beams,  one  with  a  radius  of  1  cm,  a 
peak  current  of  8  kA  and  a  current  rise  time  of  20  nsec  (NRL  Experiment)  and 
the  other  with  a  radius  of  1  cm,  a  peak  current  of  50  kA  and  a  current  rise 
time  of  60  nsec  (Sandia  Experiment).  These  calculations  can  be  made  readily 
for  a  short  time  (1  nsec),  which  is  sufficient  enough  to  provide  the  needed 
information  on  the  conductivity  generated  by  the  beam. 

VII. I  THE  RESIDUAL  CHANNEL  CONDUCTIVITY 

The  conductivity  of  the  channel  can  be  expressed  by 

8 

a  =  2.5  x  10  N  /V  (21) 

e  m 

where  is  the  momentum  transfer  collision  frequency  of  electrons  in  NH3« 
The  cross  section  for  the  electron  momentum  transfer  in  NH^  has  been  measured 
by  Pack,  et  al1 1  and  is  shown  in  Figure  2.  In  this  figure,  the  electron 
momentum  transfer  cross  section1^  in  Nj  is  also  shown  for  comparison. 

The  momentum  transfer  collision  frequencies  in  40  Tbrr  of  NH3  with 
electron  tenperatures  of  0.1  and  0.2  eV  are  ~  2.7  x  1011  sec”1  and  4  x  1011 


* 


m 


•  • 


L2 


sec  ,  respectively,  corresponding  to  the  experiments  of  NRL  and  Sandia.  Thus 


the  residual  conductivities  in  these  experiments  are  10”®  sec”1  (URL)  and 

1+5  - 1 

10  sec  (Sandia).  TSie  non-existence  of  conductivity  at  the  NRL  experiment 

has  been  borne7  out  by  the  measurements  at  NRL.  It  should  be  noted  ( See 

Figure  2)  that  when  the  plasma  electron  energy  increases  (when  the  electron 

!  beam  enters  the  channel)  the  momentum  transfer  cross  section  becomes  smaller 

and  hence  would  the  collision  frequency.  The  product  of  ov  varies  from  3.4  x 
—7  —a 

10  to  6  x  10  for  electron  energy  of  0.1  eV  to  1.0  eV  which  inplies  a 

®  lowering  of  v  by  a  factor  of  ~  6.  Even  if  this  factor  could  become  an  order 

m 

of  magnitude  the  residual  conductivity  in  the  Sandia  experiment  would  still 

6—1  Q  _ i 

be  ~  10  sec  ,  which  is  much  below  10  sec  . 

i 

VII. II  ELECTRON  BEAM  GENERATED  CONDUCTIVITY 

A  simple  hand  calculation  for  the  beam  ionization  in  ammonia  for  the  NRL 
'  and  Sandia  Experiments  is  shown  in  Figure  3  indicating  that  the  beam  generates 

much  more  electrons,  and  in  a  short  time,  conpared  to  the  residual  electron 
density.  These  calculations  are  based  on  the  assumption  that  E_  is  constant 

| 

for  times  shorter  than  1  nsec,  and  that  I  »  Iv.  With  this  the  electric  field 

n  o 

for  the  Sandia  Experiment  is  twice  that  of  NRL. 

In  Figure  4  we  show  the  electron  procuction  rates  for  direct  and 

I 

avalanche  processes  in  both  NRL  and  Sandia  Experiments,  based  on  a  more 
sophisticated  model1'*.  This  model  solves  the  rate  equations  for  ionization 
coupled  to  a  single  circuit  equation  for  Ez.  These  results  show  that  the  role 
of  avalanche  ionization  is  more  discernible  in  the  Sandia  Experiment  than  that 
of  NRL.  Also  if  the  channel  is  assumed  to  be  rarefied  (due  to  heating)  by  a 
,  factor  of  6,  the  avalanche  ionization  in  Sandia  Experiment  plays  an  inportant 

I 

role  concared  to  the  direct  beam  ionization. 


TIME  (sec) 

Figure  3  The  electron  density  build  up  in  40  Ttorr  of  ammonia  due  to  beam 
interaction.  The  avalanche  contribution  is  shown  separately  for  two 
experiments  (NRL  and  Sandia). 
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Figure  4  The  electron  density  production  rate  in  40  Torr  of  ammonia.  Direct 
beam  collisional  ionization  and  avalanche  contributions  are 
indicated.  The  density  reduction  by  a  factor  of  6  is  also  shown  for 
the  Jandia  Experiment. 


In  conclusion  one  could  say  that  in  both  experiments  the  beam  does  its 
own  thing  and  that  the  small  or  non  existent  residual  conductivity  in  the 
laser  heated  channel  adds  very  little  to  the  beam  generated  conductivity. 
However,  the  neutral  density  channel  left  behind  by  the  laser  will  modify  the 
conductivity  generated  by  the  beam. 
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